Evolution of Fluctuations during Graceful Exit in String Cosmology 
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We study the evolution of fluctuations in a universe dominated by a scalar field coupled to the 
Gauss-Bonnet term. During the graceful exit, we found non-negligible enhancements of both cur- 
vature perturbation and gravitational wave in the long wavelength limit, and we also found a short 
wavelength instability for steep background superinflation just after the completion of the graceful 
exit. This result for one possible graceful exit mechanism would provide a significant implication on 
the primordial spectrum from the string cosmology. 
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1. Introduction 
Recent developments of non-perturbative string theory 
brought a new viewpoint in cosmology. Pre-big bang 
scenario based on the scale factor duality was invented 
to predict the universe before the big bang. This scenario 
is expected to be confirmed or excluded by the observa- 
tions, where spectrum generated in the very early uni- 
verse serves as a key test. In the simple dilaton-metric 
pre-big bang model the superinflation driven by the 
kinetic energy of the dilaton is known to generate a spec- 
trum with n ~ 4 which differs significantly from the 
observationally supported n ~ 1 spectrum. It is argued 
that other field contents, e.g. pseudo-scalar axion, is re- 
quired to be included One of the most exciting pos- 
sibility pointed out in string cosmology is that the grav- 
itational waves from the stringy era might be detectable 
with the next generation gravitational wave interferome- 
ters, such as the Laser Interferometric Gravitational Ob- 
servatory (LIGO). The dilaton-driven inflation of string 
cosmology allows the existence of the stochastic back- 
ground gravitational wave which passes the constraint 
from the Cosmic Background Explorer (COBE), and at 
the same time in the detectable region of the next gen- 
eration interferometers If actually detected, this 
is supposed to be a strong evidence that the very early 
universe is described by the string theory. 

All of these discussions, however, resort to one assump- 
tion: the constancy of the long wavelength perturbations 
during the graceful exit. Since the post-big bang cos- 
mological solution corresponding to our universe is sepa- 
rated by a curvature singularity from the pre-big bang so- 
lution, it is not clear whether the fluctuations generated 
in the pre-big bang scenario is viable for the structure 
formation. 

This notorious graceful exit problem of the string cos- 
mology H is expected to be solved by non-perturbative 
quantum effects. Although the concrete mechanism is 
currently not known, there are several examples of the 
graceful exit accomplished by the quantum effects. One 



of them is a treatment using the Wheeler-de Witt equa- 
tion, and it is shown that the transition amplitude from 
the pre-big bang universe to the Friedmann universe is 
not zero Q . There is another example by Rey || , which 
is a two-dimensional dilaton-gravity model with a quan- 
tum back reaction. In this model the quantum effect 
is calculable through the conformal anomaly. There are 
several four dimensional cosmological models with higher 
curvature terms, which realize a non-singular transition 
P,M. Among them the model proposed by Antoniadis, 
et. al. |ll]] is considered as a four dimensional counter- 
part of Rey's model. This is based on the one-loop effec- 
tive action of orbifold-compactified heterotic string Jl2| , 
and consists of dilaton and modulus fields besides gravity. 
The solution in Einstein frame has an accelerated expan- 
sion phase driven by the kinetic energy of the modulus. 
This solution makes a smooth transition to a Friedmann- 
like expansion by virtue of the modulus field coupled to 
the Gauss-Bonnet curvature invariant. This is not lit- 
erally a pre-big bang model since the modulus rather 
than the dilaton drives the accelerated expansion, and 
since this acceleration is realized in Einstein frame, not 
in dilaton frame. However, this model realizes a similar 
situation to the pre-big bang in that the universe begins 
from a pole-like inflation and ends in the Friedmann uni- 
verse. Naively one can expect that the transition phase 
does not affect the superhorizon scale fluctuation, as is 
assumed in the pre-big bang scenario. 

The purpose of this letter is to study the evolution 
of the perturbations in the model of Antoniadis et.al., 
and see whether the constancy of the long wavelength 
perturbation really holds during the graceful exit. We 
hope the outcome of our study will have implication to 
the study of pre-big bang if it has a successful graceful 
exit mechanism similar to this. 

2. Graceful Exit via Gauss-Bonnet term 
The model we treat in this letter is the universe with a 
scalar field coupled to the Gauss-Bonnet invariant, which 
is essentially the same as the aforementioned one-loop 
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model |Tl|] except that the dilaton field is neglected. It 
should be noted that in several string models the coupling 
of the scalar field to the gravity we use here is shown to be 
universal beyond one-loop if there exists supersymmetry 
p3| . Thus, assuming the domination of one scalar field 
coupled to Rqb, the situation we are presenting here is 
not limited to one special model. The action of our model 
is @ 
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(1) 



where R 2 GB = R^xR/""** - $R liV B! iV + R 2 . Here, A 
is a coupling constant determined by the 4-dimensional 
trace anomaly of TV = 2 sectors, and can take both signs 
according to the relative numbers of supermultiplets in 
the theory. It is shown in pT| , p^ , ^5[ that non-singular 
cosmological solutions appear when A > 0. Since the 
absolute value of A is absorbed in the redefinition of the 
time scale, in what follows we put A = 2 for simplicity. 
For the one- loop string effective action the function £ (<p) 
is an even function of the modulus field p written in terms 
of Dedekind rj function as 



-ln[2eV( 



= -ln2-^+^--4]Tln(l 
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The existence of the non-singular cosmological solutions 
is shown by Rizos and Tamvakis E^] for general £(tp), as 
long as £ grows like or faster than p 2 as tp — > ±oo. In 
this letter we mainly concentrate on the coupling of the 
form (|J). From the action with the flat Friedmann- 
Robertson- Walker metric 



ds 2 



-dt 2 + a 2 8ijdx l dx\ 
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where a is the scale factor of the universe, one obtains 
the following equations for the background: 



if 2 =6H 2 (1-H0, 

(2H + 5H 2 )(l - Hi) + H 2 (l - £) 

Cp + 3H<p + 3(H + H 2 )H 2 ^ V = 0. 



(4) 
(5) 
(6) 



H = a/a is the Hubble parameter and the overdot de- 
notes the differentiation with respect to the cosmic time 
t. These background solutions are well studied in 0,0, 
and the numerical analysis of this model for £ = ip n is 
found in [Tq| . An example of non-singular background 
solutions (the Hubble parameter and the scalar field) 
is shown in Fig. 1. The general behavior of the non- 
singular solutions which behave like a Friedmann uni- 
verse in late time is that it starts from a Minkowski space, 
subsequently superinflates as H ~ (— i)~ 2 , and makes a 
transition at ip ~ and then goes into a Friedmann-likc 
universe (H ~ 1 /3t) . Since H and tp keep their signs 
through the evolution and the system is invariant under 



ip <-> —(p, we hereafter set H > and ip > 0. The asymp- 
totic behaviors of the modulus field are <p ~ — 51n|t| 
as i ^ — oo and p ~ y/2 /3 In t as f ^ oo. The scale 
factor approaches a constant a ~ a p in the asymptotic 
past, and grows asa~ t 1 / 3 in the asymptotic future. Of 
course, the asymptotic solutions may differ from this if we 
consider more realistic model with other matter contents 
which dominate the system in the asymptotic regions. 
If = p 2 , the past asymptotic solution becomes de 
Sitter-like, though the behavior near the transition is sim- 
ilar. 




FIG. 1. Habble parameter and scalar field in a non-singular 
solution. The accelerated expansion in the past region is 
driven by the kinetic energy of the modulus. The universe 
makes a transition to Friedmann-like solution around tp ~ 0, 
where the modulus slows down. 

3. Scalar perturbation 
To calculate the fluctuation we expand the metric and 
the modulus field up to linear order of perturbation: 

ds 2 = -(1 + 2(f)dt 2 + 2a{V\i - S i )dtdx i 

+a 2 { (1 - 2V>hy + 2£|y + 2F A] + hij } dx l dx ] , (7) 
- Sp. (8) 



p = 



We will omit the superscript of ip(°> in the following. Due 
to the symmetry of the background geometry we can de- 
compose the metric perturbations into scalar (0, ip, V, 
E), vector [Si, Fi) and tensor [hij) parts. All of these 
perturbations are assumed to be small compared to the 
background of order unity. The equations for the scalar 
perturbation becomes particularly simple if we use the 
uniform field gauge (Sp = 0, E = 0): 
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3H(3a-l)if)-6H 2 (l-a)<p 

V 2 V 2 
-2a— i; + H(3a-l) — V = 0, (9) 

(3a-l)H(f> + 2aiP = 0, (10) 
aa {<f> + (5 - 3T)iP} + (a 2 aV)' = 0, (11) 

which are time-time, time-space, traceless space-space 
part of the Einstein equations for perturbations, a and 
r are functions of the background variables, defined as 
a = 1 - Hi and T = -2H/3H 2 . The physical mean- 
ing of a is (proportional to) the fraction of the modulus 
kinetic energy to the geometrical kinetic energy, which 
can be seen from (^). T is the effective adiabatic in- 
dex F — 1 = — G l i/3G°o, where G% is the background 
Einstein tensor. They behave as a ~ t 2 , T ~ t in the 
past asymptotic region and a ~ 1, T ~ 2 in the future 
asymptotic region. In the uniform field gauge, ip is the 
curvature perturbation in terms of which the equation 
becomes 1T7I1 



A ■ V 2 

$ + + B-^i/j = 0, 



where 



A = 



a 3 a(5a 2 - 2a + 1) 
(3a - l) 2 ! 



B= (a-l) 2 (3r-2) l 
5a 2 - 2a + 1 



(12) 

(13) 
(14) 



In the long wavelength limit, the spatial derivative term 
in Eq. ( |l2| ) can be neglected and the equation is solved as 



f dt 

iP = C s (x)+D s {x) jf -. 



(15) 



Cs(x) and D$(x) are functions of spatial coordinates. 
What is expected from the analogy with the ordinary in- 
flationary universe is that, as the universe expands, the 
decrease of a -3 factor in 1/A keeps the integrand of Eq. 
( |l5| ) sufficiently small so that the long wavelength curva- 
ture perturbation stays constant. However, this is not the 
case for our model. The decrease of a (i.e. the conversion 
of the modulus kinetic energy to the universe expansion) 
overcomes the a -3 factor, which makes the integral term 
of ( Jl5| ) no longer negligible. This can be already seen by 
examining the asymptotic behaviors. In the asymptotic 
past, since A ~ t 2 , ip ~ l/(— t) and in the asymptotic 
future, A ~ t, hence if) ~ hit. In our model the variable 
which is constant in the long wavelength limit in these 
asymptotic regions is <f), the Newton potential in the uni- 
form field gauge. This cf>, however, is not always constant 
but makes a leap-like growth around the Hubble peak. 
This leap is a result of the sharp negative friction A/A 
around the Hubble peak. During the transition from the 
kinetic-driven inflation to the Friedmann universe, the 
violent movement of the background tosses up the per- 
turbation of all scales. Fig. 2 shows the behaviors of ip 
and (f> in the long wavelength limit. 



For perturbations of non-zero wavenumber k, the 
growth rate becomes more sensitive to the background 
behavior. As B is always negative for mild superinfla- 
tion, the enhancement of the short wavelength pertur- 
bation is suppressed relatively to the long wavelength. 
For steep superinflation there appears around the Hubble 
peak a period where B becomes positive, which causes 
an instability for short wavelength perturbations. The 
time evolutions of function B for different background 
initial conditions are plotted in Fig. 3. We character- 
ized the background behavior by the peak value of the 
Hubble parameter, i/max- While the growth of pertur- 
bation caused by A/ A is independent of the wavenumber 
k, the effect of B depends on the scales of the perturba- 
tions. The transfer function of the spectrum during the 
graceful exit will thus depend largely on the background 
behavior, i.e. the initial conditions. 

We have to point out the possibility that in realistic 
string theory higher curvature contribution might sup- 
press this small scale instability. Since this instability 
occurs for large R ~ H 2 , our assumption of R% B dom- 
inance might not be valid. Thus, we cannot conclude 
from this analysis alone that this instability is inherent 
in the string theory. 



FIG. 2. Time evolution of long wavelength perturbations 
(a) ip and (b) <f>. The background is the same as Fig. 1 
(JJmax = 0.5). The normalization of the amplitude is set as 
d> = 1 at t = 0. 
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4- Rotation and Gravitational Wave 
From the vector part of the perturbative equation we 
can see the behavior of the vorticity. The equation can 
be solved to give Si + aFi = 0. This is a consequence of 
the absence of the rotational sources. 

The tensor part of the perturbative equation is 



3H 



h 



V 2 

— (5-3T)/i 
ar 



0. 



(16) 



where h is the Fourier component of either plus- or cross- 
mode gravitational wave amplitude. In the long wave- 
length limit the equation can be solve as 



r* dt 

h = C T (x) + D T (x) / 



(17) 



and the second term (integral) represents the growth of 
the perturbation which is very similar to ifj of the scalar 
mode. The short wavelength gravitational wave is quite 
different from the scalar mode perturbation. Since 5 — 3r 
is always positive before the Hubble peak, the gravita- 
tional wave is exponentially enhanced from the asymp- 
totic past. To avoid the breakdown of the background 
solution, initial perturbation has to be fine-tuned to be 
small Hff9|. 



H_max=10 




3. Time evolutions of B for different initial conditions. 
For mild superinflation B remains negative during the tran- 
sition. As the superinflation becomes steep, there appears a 
period where B is positive, leading to an unstable mode for 
the short wavelength perturbations. The figure shows three 
examples, specified by the maximal values of the Hubble pa- 
rameter, i/max = 0.1, 1, 10 (-ffmax is larger for steeper su- 
perinflation). For comparison t = is defined to be the time 
when H = i/max, i.e. H = 0. 



supposed to be more reliable choice of the gauge in order 
to avoid the breakdown of the perturbations. It might be 
suspected that the growth of the fluctuations we found 
here might be of this kind, which can be "gauged down." 
However, in our model the choice of off-diagonal gauge 
or longitudinal gauge will lead to more violent growth of 
the perturbation, and as can be seen from the constancy 
of <j) in the asymptotic regions, the uniform field gauge we 
used here is supposed to be one of the most gentle gauge 
choice. While the cancellation of the growth may occur 
in a very accidental case, the growth we found here is a 
generic result and possibility of gauging down is unlikely. 

Although our study is limited to only one simple 
model, there are reasons to believe this is a quite general 
phenomena that may arise in string-cosmological mod- 
els with a successful graceful exit. One reason is that 
the coupling we used here, Eq. (|^), appears in many 
string low-energy effective actions as a loop correction. 
Moreover, recent studies based on duality show that the 
function of the form (|J) appears as non-perturbative R 2 
couplings suggesting that this coupling form is ex- 
act beyond one-loop. The extension of our analysis to 
multiple fields, 



S 



A 

To 



i{Vi)R 2 GB 



(18) 



leads to essentially the same equations of motion for per- 
turbations ^0|. The enhancement of the fluctuations 
for scalar and tensor modes is a consequence of non- 
negligible l/a 3 a, which means the rapid conversion of 
the scalar field kinetic energy into the kinetic energy of 
the universe expansion. Since this situation is the same 
for any cosmological model with a transition from a su- 
perinflation to Friedmann universe, we expect the ex- 
istence of the same enhancement of the perturbations 
even for more realistic models with dilaton or other fields, 
as long as the graceful exit is accomplished through the 
same mechanism. The model-independent approach by 
Brustein et.al. ]2l| is of interest from this viewpoint. 

We are grateful to Jai-chan Hwang and Hyerin Noh 
for the discussions in 19th Texas Symposium in Paris. 
We also appreciate the useful discussions with Atsushi 
Taruya. Part of this work by J.S. is supported by the 
Grant-in- Aid for Scientific Research No. 10740118. 



5. Discussions 
We have shown that during the transition there exist the 
enhancements of perturbations in a cosmological model 
dominated by a scalar field coupled to the Gauss-Bonnet 
term. In the case of the dilaton-driven pre-big bang 
model, there is a growth of perturbations as the universe 
evolves towards the curvature singularity. While this 
growth is present in the longitudinal gauge, it is strongly 
suppressed if one uses the off-diagonal gauge Q , which is 
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